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Abstract: This paper presents a novel method based on a fiducial marker 
for correction of motion artifacts in 3D, in vivo, optical coherence 
tomography (OCT) scans of human skin and skin scars. The efficacy of this 
method was compared against a standard cross-correlation intensity-based 
registration method. With a fiducial marker adhered to the skin, OCT scans 
were acquired using two imaging protocols: direct imaging from air into 
tissue; and imaging through ultrasound gel into tissue, which minimized the 
refractive index mismatch at the tissue surface. The registration methods 
were assessed with data from both imaging protocols and showed reduced 
distortion of skin features due to motion. The fiducial-based method was 
found to be more accurate and robust, with an average RMS error below 20 
|j,m and success rate above 90%. In contrast, the intensity-based method had 
an average RMS error ranging from 36 to 45 |j,m, and a success rate from 
50% to 86%. The intensity-based algorithm was found to be particularly 
confounded by corrugations in the skin. By contrast, tissue features did not 
affect the fiducial-based method, as the motion correction was based on 
dehneation of the flat fiducial marker. The average computation time for the 
fiducial-based algorithm was approximately 21 times less than for the 
intensity-based algorithm. 
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OCIS codes: (170.4500) Optical coherence tomography; (100.6950) Tomographic image 
processing; (100.2000) Digital image processing. 
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1. Introduction 

Optical coherence tomography (OCT) [1] provides a means to assess skin pathologies non- 
invasively. Its application in dermatological research and real-time, in situ examination of 
skin and skin diseases in clinical scenarios has been explored in recent years [2,3]. Skin 
structures identified using OCT, such as the stratum comeum, epidermis, dermo-epidermal 
junction, blood vessels, hair follicles, sweat ducts, and sebaceous glands, have been found to 
show good correlation with histology [1,4]. OCT has been utilized to detect skin cancer [5], 
monitor skin changes by UV radiation [6] or aging [7], examine skin inflammatory diseases 
[3] and vasculature dilation [3,8], as well as the efficacy of laser treatment for vascular 
malformations [4]. 

Common practice in dermatological application of OCT is the use of a handheld imaging 
probe for convenient access to different body sites [9-13]. At the same time, there has been 
growing interest in utilizing three-dimensional OCT (3D-0CT) to aid in visualization and 
improve clinical assessment. Such 3D-OCT scans are commonly constructed from a series of 
2D vertical cross-sectional images (B-scans) acquired by raster scanning across the skin. 
Typically, such scans are subject to significant motion artifacts, due to cardiac and respiratory 
motion, and to inadvertent shaking of the handheld probe. 

Our experience has been that motion artifacts are minimally present within each individual 
2D B-scan but between consecutively acquired B-scans, relative offsets and rotations distort 
the appearance and size of anatomical features, such as the skin's surface topology and 
subsurface blood vessels. The presence of motion artifacts as a source of 3D-OCT image 
degradation has been reported in previous publications [9,14,15] involving in vivo human 
skin. In several publications [5,16,17], it has also been indicated that a portion of the clinical 
images of skin have been omitted from analysis due to motion artifacts. 

Previous efforts in compensating motion artifacts in OCT images have been reported in 
the fields of ophthalmology [18-21], endoscopy [22,23], and cardiology [24]. Hardware 
techniques to reduce motion artifacts during in vivo imaging include: decreasing image 
acquisition time by using a higher scan rate (such as megahertz OCT [25]); avoiding dense 
sampling; or restricting scanning to a small area of the sample. 

An alternative compensation approach adopted in other imaging modalities is post- 
processing using image registration [26], which is the process of aligning two or more images 
into spatial correspondence. Manual image registration is highly time-consuming. For a 3D- 
OCT image, this can involve the manual correction of hundreds or thousands of B-scans. 
Fully automated image registration methods are necessary in order to incorporate such 
correction into a clinical workflow. Automation also removes the subjectivity inherent in 
manual correction. 
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There are two main categories of image registration methods: intensity-based and feature- 
based methods. Intensity-based methods ahgn images based on their pixel intensity [27]. Such 
methods involve applying a geometrical transformation to one image (referred to as the 
floating image), and calculating the pixel-by-pixel similarity to the other (referred to as the 
fixed image). The geometrical transformation is iteratively modified to optimize the similarity 
of the floating and fixed images. This similarity is quantified with a 'similarity measure' [27], 
a single value calculated from the intensity values of pairs of corresponding pixels from the 
two images. Choice of similarity measure, resampling method during the application of the 
geometrical transformation, and optimization strategy constitute the main differences between 
various intensity-based methods. As an example, Zawadzki et al. [20] utilized cross- 
correlation in an intensity-based registration to ahgn consecutive B-scans in a 3D-OCT image 
of the eye. Kraus et al. [28], in contrast, utilized the sum of the squared L2 norm as a 
similarity measure to correct eye motion in multiple, orthogonally acquired volumetric scans. 

By contrast, in feature-based methods, the geometrical transformation between two images 
is calculated by extracting specific tissue structural features. These features are used to 
establish the correspondence between the two images. Examples of such features include 
lines, curvature extrema, or contours extracted from each B-scEui [29,30]. Antony et al. [21] 
described a surface segmentation-based registration method to correct for distortions seen 
along the fast- and slow-scanning axes in 3D-0CT retinal scans. 

Feature-based methods may also include the use of a fiducial marker, an object with a well 
characterized size and shape, to provide features with which to register images. In this paper, 
we present a fiducial marker specifically designed for OCT imaging of skin. The fiducial is a 
small, flat metal square with a circular hole through which the skin may be imaged. By 
affixing the fiducial marker to the skin, the flat surface of the fiducial may be automatically 
identified and used to correct motion between adjacent B-scans. 

In our OCT imaging of skin, two different imaging protocols are used based on whether or 
not a refractive index (Rl)-matching medium is applied. Use of a Rl-matching medium has 
been shown to reduce the impact of a range of morphological and intensity artifacts common 
in OCT imaging of skin [31]. Application of such media to the skin has also been shown to 
modify the measured thickness of the epidermis. Thus, the choice of protocol is dependent 
upon the specific clinical application. 

Whilst each imaging protocol is subject to motion artifact, it manifests in different ways. 
As described in [31], when matching the RI of skin, it is typical to apply a substance such as 
ultrasound gel and then hghtly compress the skin under a glass cover slip, eliminating air from 
the imaging surface. Because of this contact with the skin, motion artifacts are predominantly 
perpendicular to the axial direction of the light beam. When imaging directly from air into the 
skin, no glass cover slip is necessary and motion artifact between B-scans tends to be 
primarily in the axial direction. In both situations, slight rotations occur between adjacent B- 
scans. 

The goal of the present study was to develop a fully automated motion correction scheme 
based on the novel use of a fiducial marker, and compare this method against a standard 
intensity-based registration approach. We have assessed these methods both with and without 
Rl-matching medium. Comparison of the methods was carried out in terms of registration 
accuracy, robustness and speed, and assessed against a gold-standard manual image 
registration. This is the first pubhshed assessment of such methods for in vivo 3D-OCT scans 
of human skin. 

2. Materials and methods 

2.1. Data 

The image registration algorithms were assessed on both normal and scarred skin, 
representing two distinctly different tissue sub-types in terms of surface texture and 
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subsurface structure. The wound healing process in skin resuhs in a broad spectrum of scar 
types, ranging from normal fine-line scars to various abnormal and problem types that require 
treatment, such as hypertrophic and keloid scars, and scar contractures [32]. The surface 
texture of scars ranges from flat and smooth surfaces to raised, irregular surfaces with large 
corrugations. Unlike normal skin, such scars seldom have a regular pattern of surface furrows 
and may be highly perfused by blood vessels [33] causing them to appear red or inflamed. 

Four subjects with scars and two healthy subjects (four males, two females, mean age: 32 
years) were enrolled with prior informed consent. Hair on the skin/scar area to be imaged was 
trimmed using an electric shaver prior to scanning. Each subject was scanned in 1-3 locations, 
giving a total of 7 scar and 7 normal skin locations across all subjects. For each of the 14 
locations, two 3D-0CT scans were performed, giving a total of 28 data sets. The two scans 
acquired at each location are illustrated in Figs. 1(c) and 1(d): one with the skin exposed to air 
and one with a thin layer of ultrasound gel (350 ± 100 (xm thick) sandwiched between the 
glass cover slip and the skin/scar. Each 3D-0CT scan was 6.5 x 6.5 x 3 mm (x x j x z) in size 
and comprised 1088 x 1088 x 512 pixels. We define the x andy axes to represent the lateral 
dimensions of the scan, and z to represent the axial depth dimension. Raster scanning was 
performed with fast-axis B-scans orientated parallel to the xz plane, and slow-axis B-scans 
parallel to the yz plane. The acquisition time per 3D scan was 2 minutes and 20 seconds. 




Fig. 1. Handheld SS-OCT probe and sample spacer setups for in vivo imaging of skin, (a) 
Photograph of the imaging probe setup showing interchangeable sample spacers, (b) 
Photograph of the handheld probe in use. (c) Schematic diagram of imaging protocol in air. (d) 
Schematic diagram of imaging protocol in Rl-matching medium (ultrasound gel). 



2.2. OCT system 

A fiber-based swept-source polarization-sensitive OCT system (PSOCT-1300, Thorlabs, 
USA) comprising a broadband swept laser source and a Michelson interferometer with 
balanced detection was used in this study. The laser source has a central wavelength of 1325 
nm and the average optical power of the probing beam incident on the skin was 3.2 mW. Data 
was acquired using in-built software speckle reduction, which reduced the system's axial 
resolution. The resulting FWHM axial resolution was measured to be 20 fim (in free space) 
and the lateral resolution was measured to be 16 (xm. The working distance and numerical 
aperture of the objective lens are 25.1 mm and 0.056, respectively. The system has an axial 
scan rate of 16 kHz. 

The sample arm was configured as a handheld imaging probe, whereby one of two 
adjustable sample spacers could be affixed to the probe (Fig. 1) to maintain a constant 
distance between the objective lens and the skin. Spacer 1 consists of a 25 mm 0 hole to allow 
scanning of the skin in air (i.e., no RJ-matching) whilst the probe is held firmly in contact with 
the skin. Spacer 2 has a smaller hole (18 mm 0) and a grooved base in order to attach a glass 
cover slip (22 x 40 x 0.13 mm), which was used as a window for imaging the skin through a 
layer of Rl-matching medium (ultrasound gel). 
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2.3. Image acquisition 

A small metal fiducial marker (inset of Fig. 1(a)) was adhered securely onto the skin using 
double-sided adhesive tape. The marker is a brass shim (1 cm square, thickness: 170 (xm) with 
a 5 mm diameter imaging hole. The marker was used to perform the gold-standard manual 
image registration and facilitate the fully automated feature-based registration algorithm. 
During image acquisition, the imaging probe was placed in contact with the subject's skin and 
a 3D scan was acquired across an area covering the entire hole of the fiducial marker (Fig. 1). 
To reduce the intensity of specular back reflections to the probe, the axis of the scanning beam 
was tilted at a small angle (~1°) to the surface normal of the marker and glass cover slip. The 
OCT beam focus was set at -200 |a,m below the skin surface for optimal subsurface imaging. 

2.4. Image registration methods 

The motion artifact between adjacent fast-axis B-scans was modeled with a rigid 
transformation consisting of a rotation and a 2D translation. Transformation of each B-scan to 
the first reference B-scan was subsequently computed by concatenating adjacent 
transformations, allowing reconstruction of the motion-corrected 3D data set. The two 
registration approaches we have adopted to estimate the adjacent geometrical transformations 
are detailed in the following sections and illustrated in Fig. 2. 

2.4.1. Feature-based registration using the fiducial marker 

The intersection of the fiducial marker with a single B-scan imaging plane appears as a 
strongly reflecting horizontal line with a large break in the center, corresponding to the 
surface of the fiducial and the circular imaging area, respectively. Thus, the feature-based 
method involved two steps: automatic segmentation of the fiducial marker, which appears as 
two line segments in each B-scan; and registration of these line segments between adjacent B- 
scans. 



Acquisition of 
3D data set 



Manual segmentation 
of tissue region 

(En face average image ) 

y Discarded 



Registration accuracy assessment 




Automated 
feature-based registration 

(1) Segmentation 



Automated 
intensity-based registration 




Image A Image B 

(2) Registration (smoothed) (smootRed) 

Image B , ' 




Adjacent transfonnatlon 
parameters (0, tx and t^) 



Computation of the transformation 
parameters for each B-scan 
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Transformation of 
the 3D data set 
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— Fitted model • Before registration 
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Fig. 2. Flow diagram for tlie comparison of feature-based and intensity-based registration 
methods. 



#167671 - $15.00 USD Received 2 May 2012; revised 21 Jun 2012; accepted 26 Jun 2012; published 29 Jun 2012 
(C) 2012 OSA 1 August 2012 / Vol. 3, No. 8 / BIOMEDICAL OPTICS EXPRESS 1779 



Line segmentation used the Hough transform [34], which is a standard image processing 
algorithm widely used in other imaging modalities [35,36]. Results are shown in Fig. 2, as a 
magenta line overlay on a B-scan (far left image in the middle row). Prior to applying the 
Hough transform, 20 pixels at the start of each A-scan were cropped to eliminate the 
confounding effects caused by autocorrelation noise. As the metal fiducial gave a strong 
reflection, a threshold (Ti) was applied to mask out other features (e.g., skin). The Hough 
transform was then performed on the left l/8th and right l/8th of the image. The central 
portions of the image typically contained skin and not fiducial, and were not included in the 
computation. 

During execution of the Hough transform, an accumulator space (i.e., a 2D histogram) was 
created and quantized into an array of (p, 0) bins, in which p is the perpendicular distance 
from the line to the origin, and 0 is the angle that the perpendicular line makes with the x-axis, 
as defined in [34]. The global maximum of the histogram was identified, and all connected 
histogram bins with a high count (>95% of the global maximum) were tagged. The algorithm 
was then used to calculate the average (p, 0) values of the tagged histogram bins, giving the 
parameters of the line corresponding to the surface of the fiducial marker. 

The inner edges (i.e., end points) of each line segment (left and right) of the marker were 
then identified by tracing along the line until the image intensity fell below an empirically 
chosen threshold value, T2. At conclusion of the calculation, two line segments had been 
exfracted, corresponding to the intersection of the fiducial marker with the B-scan on each 
side of the image. 

To register adjacent B-scans, the rotation angle, 6ab between them was determined from 



where A and B refer to the fixed and floating images, respectively, and tha and thb are the 
slopes of the detected lines in each pair of adjacent B-scans. The horizontal and vertical 
ttanslations between B-scans were determined from Eq. (2): 



where and are the horizontal and vertical franslations, respectively, (xa, za) and (xb, zb) are, 
respectively, the coordinates of the cental point of the line segments detected in images A and 
B. 

A range of threshold values (Ti and T2) for segmentation were tested and optimal values 
were chosen empirically to reduce segmentation error. The correct choice of thresholds may 
affect the success of a segmentation process (and, thus, the feature-based regisfration process) 
by affecting the accuracy of detection of the fiducial marker's surface and edges. Robustness 
of the algorithm over a range of values for Ti and T2 will be discussed in Section 4. 

2.4.2 Intensity-based registration 

For each pair of adjacent B-scans, the intensity-based method searches for an optimal 
geomefrical mapping to transform the floating image onto the flxed image. We use the 
correlation coefficient as the similarity measure between the images [26,27]. The optimal 
geomefrical transformation will maximize the correlation coefficient between the images, 
defined as 



Eq. (1): 




(1) 




CC = 



Z..n^ (A(x,)-A)(B'(x,)-B) 



(3) 
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where A and B denote the fixed and floating images, respectively; is the pixel position in 
image A, A{xa) and B\xa) are the pixel intensity values of these images at position x/, A is the 
mean pixel intensity value in image A and B is the mean of L , where B^ is image B 

transformed to the coordinate space of image A, and CIj^b is the overlap domain of images A 
and fl^ for a given transformation estimate T. 

An iterative search for the optimum transformation at sub-pixel resolution was performed 
using the pattern search optimization method [37]. The search was initialized with the identity 
transformation. At each iteration, Harming- windowed sine interpolation (half width: 18 (xm) 
was applied to resample the floating image to the coordinate space of the fixed image. A 
common source of bias in registration algorithms is the failure to account for the varying 
degree of overlap of the images with each geometrical transformation, which impacts the 
numerical value of the similarity measure. A fixed number of overlap pixels (within a 
rectangular window) between both images was used for the correlation coefficient calculation 
in order to avoid such bias. 

To avoid the influence of the fiducial marker on the intensity-based registration, the skin 
region of each data set was delineated before registration. Cropping of the marker allowed 
simulation of the situation in which a fully automated intensity-based registration is 
performed on 3D data sets without the presence of the marker. This delineation was 
performed manually to avoid bias from segmentation errors. 

To reduce the confounding effects of speckle noise on the intensity-based registration 
algorithm, each fast-axis B-scan was low-pass filtered prior to registration using a 2D 
anisotropic Gaussian kernel. Several smoothing factors (defined as the ratio of FWHM of the 
Gaussian smoothing kernel to that of the measured system resolution) were tested and the 
registration results obtained using the optimum factor are reported. 

2.5. Technique for comparison of registration accuracy and robustness 

Figure 2 illustrates the method used to compare the registration algorithms. A gold-standard 
manual registration was performed on all 28 data sets by aligning the fiducial marker visible 
within each B-scan. The accuracy of each registration algorithm was quantified at 50 
manually selected control points distributed throughout each data set. These control points 
were positioned on the edge of the circular imaging field (the inner boundary of the marker). 
For each registered data set, we calculated the minimum RMS error (in microns) between 
each of the 50 geometrically transformed control points and the closest point on a manually 
fitted model of the inner boundary of the marker. 

To assess the robustness of each registration algorithm, we calculated the RMS error both 
before and after its application. Where the RMS error was reduced, the algorithm was deemed 
successful, whilst an increase in RMS error was recorded as a failure. To assess accuracy, the 
average RMS error was calculated for data sets for which an algorithm was deemed 
successful. 

Computation times for the registration algorithms were also recorded. The algorithms 
were implemented in MATLAB (vR2009b, Mathworks, Natick, Massachusetts) on an Intel 
Quad Core i7@3.07GHz computer. The average computation time to register a 3D volume 
consisting of approximately 750 ± 30 fast-axis B-scans using the feature-based method, 
including the time required for segmentation, was 2.4 min. This is approximately 21 times 
faster than intensity-based method, which required an average of 51 min per registration. 
Algorithm implementations were not optimized for execution speed. 

3. Results 

Table 1 shows the success rate and average RMS error obtained for both registration methods, 
grouped by imaging protocol. We will refer to the data sets collected with these imaging 
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protocols as 'air-tissue' and 'gel-tissue'. The success rate of both algorithms varies with the 
imaging protocol. 



Table 1. Registration accuracy and success rate for air-tissue and gel-tissue data sets 



Registration method 






Average RMS error, 


Success rate, 


(optimum parameter(s)) 


Data sets 


nm (±SD)* 


% 


Feature-based method 


Air-tissue (n = 


14) 


12 (±3) 


100 


(Tl: 98%; T2: 60%) 


Gel-tissue (n = 


14) 


19 (±8) 


93 


Intensity-based method 


Air-tissue (n = 


14) 


36 (±13) 


50 


(Smoothing factor: I) 


Gel-tissue (n = 


14) 


44 (±20) 


86 



*Note that only successful registrations were used in computing the accuracies (average RMS errors), n: number of 
data sets. 



The feature-based method using the fiducial marker was found to have comparable 
success rates for air-tissue data sets (100%) and gel-tissue data sets (93%) for Ti and T2 values 
of 98% and 60%, respectively. A success rate above 90% for air-tissue data sets required a Ti 
value of greater than 70%. To obtain a similar success rate for gel-tissue data sets, it was 
necessary to increase Tj to 95% in order to ehminate extraneous horizontal lines from the 
glass cover slip and artifacts caused by detector saturation. Success of the feature -based 
registration was robust imder a wide range of values for T2. A success rate of over 90% was 
achieved for 40% < T2 < 70%) under both imaging protocols. For T2 lower than 40%), the 
edges of the hole could not be determined accurately primarily because of noise in the image. 
For T2 higher than 70%), parts of the marker's surface became invisible, resulting in inaccurate 
segmentation. The average RMS error of the feature-based method was found to be 
approximately the same as the error associated with the memual selection of control points, 
which was estimated to be 15 ± 3 jxm (»2 to 3 pixels). 

The average RMS error for the intensity-based method varied between imaging protocols, 
and in both situations was inferior to the feature -based method. The intensity-based method 
gave a higher success rate for the gel-tissue protocol (86%) than for the air-tissue protocol 
(50%)). We foimd that the intensity-based method had its highest success rate when the image 
was smoothed with a smoothing factor (defined in Section 2.4.2) equal to 1. When no 
smoothing was apphed, the accuracy was adversely impacted by the presence of interpolation 
artifacts [38]. When a smoothing factor of 2 or more was applied, blurring of small features 
and edges reduced the accuracy. 

Normal and scarred skin typically have very different surface and subsurface structures. 
These structures have an influence on the success of the intensity-based method. Table 2 
provides the registration accuracy and success rate of this method organized by tissue type. 
Registration of normal skin data sets resulted in a slightly lower error (38 |iim) and higher 
success rate (79%) compared with scar data sets (45 (xm and 57%, respectively). By contrast, 
tissue features did not affect the feature-based method, as the registration was based on 
delineation of the fiducial marker. 



Table 2. Intensity-based registration accuracy and success rate for normal skin and scar 

data sets 



Data sets 


Average RMS error, |.im (±SD)* 


Success rate, % 


Normal skin (n = 14) 


38 (±11) 


79 


Scar (n = 14) 


45 (±25) 


57 



*Note that only successful registrations were used in computing the accuracies (average RMS errors). Results are 

based on the use of a smoothing factor of 1 . 



Figures 3 and 4 illustrate typical registration results for air-tissue and gel-tissue data sets, 
respectively. Figure 3(a) shows the impact of uncorrected motion artifact on an OCT 
acquisition of scarred skin. The scar surface and subsurface features such as blood vessels are 
significantly distorted. With air-tissue imaging, these distortions were primarily due to vertical 
motion. Both feature-based (Fig. 3(b)) and intensity-based (Fig. 3(c)) methods were shown to 
enhance the visualization and assessment of the scar by reducing the amoimt of motion 



#167671 - $15,00 USD Received 2 May 2012; revised 21 Jun 2012; accepted 26 Jun 2012; published 29 Jun 2012 
(C) 2012 OSA 1 August 2012 / Vol. 3, No. 8 / BIOMEDICAL OPTICS EXPRESS 1782 



artifact. The third row of Fig. 3 shows a representative slow-axis B-scan (yz -plane), the 
location of which is indicated by the vertical magenta dot-dashed line in the second row of 
Fig. 3. The bottom row shows a slow-axis B-scan of the fiducial marker (at a location 
indicated by purple dot-dashed line in the second row). Errors in the intensity-based algorithm 
(right column) manifest as a curve in this line. 




Fig. 3. Typical results for an air-tissue data set: (a) before registration; (b) after feature-based 
registration; and (c) after intensity-based registration. For each of (a), (b), and (c): 3D solid 
render side view (Row 1) and top view (Row 2); a slow-axis B-scan across the scar (Row 3); 
and fiducial marker (Row 4) at positions indicated by the magenta and purple dot-dashed lines, 
respectively. 



Motion artifacts in data sets acquired with the gel-tissue imaging protocol were observed 
to occur primarily along the x and y axes, as illustrated in Fig. 4(a). The third row of Fig. 4 
shows a slow-axis B-scan acquired along the vz -plane, and two examples of motion artifact 
are delineated by the yellow dashed boxes. Both registration methods were shown to 
effectively reduce motion artifacts in the gel-tissue data sets. With this imaging protocol, the 
glass cover slip is apparent, as shown in the fourth row in Fig. 4. The additional horizontal 
lines at depths below the fiducial marker are artifacts caused by detector saturation. This is 
inevitable because the marker exhibits a very strong back reflection that is outside the 
dynamic range of our detector. This signal, as well as its beat terms with the reflecting glass 
cover slip interfaces, is clipped by the finite dynamic range of the digitizer input, resulting in 
harmonics in the FFT reconstruction. The marker's surface, nonetheless, was found to rehably 
appear as the strongest and widest horizontal line in data sets in this study acquired by tilting 
of the glass cover slip to reduce the signal strength of the extraneous hnes. 
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Fig. 4. Typical results for a gel-tissue data set: (a) before registration; (b) after feature-based 
registration; and (c) after intensity-based registration. For each of (a), (b), and (c): 3D solid 
render side view (Row 1) and top view (Row 2); a slow-axis B-scan across the scar (Row 3) 
and fiducial marker (Row 4) at positions indicated by the magenta and purple dot-dashed lines, 
respectively. Yellow dashed boxes indicate the distortion of the scar features by motion in 
horizontal x direction. Arrows in Row 4: cyan - surface of the marker; pink - upper and lower 
interface of the glass cover slip; and orange - examples of artifacts caused by detector 
saturation. 

4. Discussion 

Motion artifact is a source of image distortion that can compromise the clinical utility of in 
vivo 3D-0CT imaging. The magnitude of the motion artifacts observed in our system was on 
the order of tens to hundreds of microns between adjacent fast-axis B-scans, with more severe 
distortions observed on images acquired from areas of the body more affected by involuntary 
movement, such as the neck, chest, and abdomen. 

This study has proposed a novel method based on a fiducial marker for correction of 
motion artifacts in 3D. This feature-based approach was compared to a standard cross- 
correlation intensity-based registration method and found to be more accurate (smaller RMS 
error) and reliable (higher success rate) than the intensity-based method. It was observed that 
the two imaging protocols assessed were susceptible to different motion artifacts. In general, 
both registration methods were effective in reducing these artifacts for both imaging 
protocols, as evidenced by the smooth and continuous skin features observed in the slow-axis 
B-scan after registration (Figs. 3 and 4). 

Potential disadvantages of the feature-based algorithm are the requirement for a fiducial 
marker to be present within the imaging field and the reduced field of view. The use of a 
fiducial marker may not be feasible in some clinical situations, such as with fresh wounds. 
Our experience is that the use of the marker did not present any complications in our clinical 
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scanning conducted during follow-up burn scar assessments on successfully grafted scar 
tissue. 

The success of the feature-based method also depends heavily on the success of the 
fiducial segmentation process. The most common errors arose in the gel-tissue imaging 
protocol, in which the glass cover slip resulted in additional horizontal lines in the field of 
view. Even though the marker in general produces a much stronger signal, in some scans the 
additional horizontal lines confounded its automatic detection. Optimal choices of threshold 
parameters Ti and are important in reducing such errors. We note that the algorithm was 
found to be robust over a wide range of values for Ti (>70% for air-tissue; >95% for gel- 
tissue data sets) and T2 (40% to 70% for both imaging protocols). Such robustness suggests 
that the algorithm may be applicable to a wide range of data sets without the need for 
adjustment of these parameters. 

The intensity-based method was found to be adversely influenced by surface furrows and 
corrugations, often misinterpreting these variations in surface topology as motion artifact and 
inappropriately attempting to smooth these features. Over an entire 3D-0CT scan, we 
observed that this effect introduced a bias in the positioning of the set of B-scans, resulting in 
substantial cumulative misplacement of the final B-scan. This resulted in a reduction in the 
measured accuracy and success rate of the intensity-based algorithm. This effect was most 
pronounced in the registration of scar data sets (Table 2), as a majority of the scars in this 
study had more significant surface corrugations than were present in normal skin. The 
distortion caused by this bias was substantiated by the deformation of the fiducial marker, 
which was simultaneously registered in the registration process. Similar 'flattening' artifacts 
have been reported in the registration of OCT images in the eye [20]. Even at higher OCT 
acquisition speeds, where motion is minimized, this bias will still be present due to the skin's 
surface corrugations. 

Intensity artifacts manifest as low intensity vertical streaks specifically located under the 
skin surface fiirrows [31] and partially account for the lower success rate of the intensity- 
based algorithm with the air-tissue imaging protocol as compared to the gel-tissue protocol 
(50% versus 86%). These artifacts are more significant in the air-tissue protocol (Figs. 5(a), 
5(c), Row 2), and are reduced after the application of ultrasound gel (Figs. 5(b), 5(d), Row 2). 
Intensity-based algorithms implicitly assume a very similar distribution of intensity values 
between adjacent B-scans, and local intensity artifacts reduce the validity of this assumption. 
Figures 5(a) and 5(c) show the skew in the round hole of the fiducial marker as an effect of 
this artifact for data sets such as these with major surface furrows oriented at a diagonal angle 
to the fast-scanning axis. Figures 5(b) and 5(d) show the co-located scans acquired with the 
gel-tissue protocol, which reduced this artifact. The shape of the fiducial hole is better 
preserved in these images. 




Fig. 5. Intensity-based image registrations on normal skin ((a) and (b)) and scarred skin ((c) 
and (d)) with an air-tissue ((a) and (c)) and gel-tissue ((b) and (d)) image protocol, (a) and (b) 
are co-located as are (c) and (d). Upper row: en face surface images. Lower row: OCT B-scans 
at location indicated by the red dotted line. An'ows: black - direction of fast-scanning axis; 
white - intensity artifacts. 
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While failure of the intensity-based algorithm sometimes resulted in subtle corruption of 
the skin features, we observed that the feature-based algorithm tended to clearly succeed or 
clearly fail. Such catastrophic failures are more easily detected by manual inspection of the 
registered data set. While all algorithms will fail in a subset of situations, we consider it 
desirable in a clinical setting for such failures to be easily detectable. 

The motion correction algorithms outlined in this paper do not account for motion along 
the slow-scanning direction (y-axis). Such motion will result in uneven spacing of the B- 
scans. The well-defined geometry of the fiducial marker may allow for extensions to account 
for such motion, by assessing the shape of the flat marker in the en face plane. For example, 
with the circular marker utilized in this paper, we note that it may be feasible to utilize a 
circular Hough transform [39] to automatically extract its correct shape. 

Computation time of the registration method has important implications for the feasibility 
of 3D-OCT in on-site clinical examination. The feature-based method was found to perform 
21 times faster than intensity-based method because it operates on a small number of points 
and is not iterative. By contrast, in the intensity-based method, a similarity measure based on 
all pixel values in each pair of adjacent fast-axis B-scans is calculated, and the geometrical 
tremsform between each pair of images is optimized iteratively. This large contrast in the 
computational complexity of the two algorithms is in agreement with comparisons of feature- 
based and intensity-based registration algorithms undertaken in other imaging modalities [28]. 

5. Conclusion 

We have described a method for motion correction of 3D-0CT scans of skin using a fiducial 
marker with feature-based image registration. The method has been quantitatively compared 
against a cross-correlation intensity-based registration method, and accuracy was assessed 
against manual motion correction. The methods were tested on a total of 28 data sets acquired 
using two imaging protocols: from air into tissue; and with ultrasound gel as an optical 
coupling layer. The feature-based method was shown to be more accurate and reliable, with 
an average RMS error of below 20 (xm and success rate of above 90%. The intensity-based 
method had an average RMS error ranging from 36 to 45 (j,m and a success rate ranging from 
50% to 86%. The performance of both algorithms was found to be influenced by the imaging 
protocol. Additionally, the intensity-based method was affected by the level of corrugation of 
the skin surface. Our results show that the use of a fiducial marker can greatly aid in motion 
correction of in vivo 3D-0CT skin scans. 

Acknowledgments 

The authors gratefully acknowledge the following colleagues who contributed to this study: 
Rodney W. Kirk, Lixin Chin, Ngie Min Ung, Pei Jun Gong and Timothy R. Hillman. The 
authors also thank Sharon Rowe at the Telsfra Bums Reconsfruction and Rehabilitation Unit, 
Royal Perth Hospital, for her help in patient management. Y. M. Liew is supported by a 
scholarship from The University of Western Ausfralia. R. A. McLaughlin is supported by 
funding from the Cancer Council Western Australia, and this study is in part supported by the 
National Breast Cancer Foundation, Austtaha. 



#167671 - $15,00 USD Received 2 May 2012; revised 21 Jun 2012; accepted 26 Jun 2012; published 29 Jun 2012 
(C) 2012 OSA 1 August 2012 / Vol. 3, No. 8 / BIOMEDICAL OPTICS EXPRESS 1786 



